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Nanofluids are considered to have great potential for heat transfer enhancement and are highly suited 
to application in practical heat transfer processes. Recently, several important studies were carried out 
to understand and explain the causes of the enhancement or control of heat transfer using nanofluids. 
The main aim upon which the present work is based is to give a comprehensive review on the research 
progress on the natural convective heat transfer characteristics of nanofluids for both single- and two- 
phase models. Both experimental and theoretical studies are reviewed for natural convection of 
nanofluids in different types of enclosures. 
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1. Introduction 


Natural convection heat transfer is an important phenomenon in 
engineering and industry with widespread application in diverse 
field, such as, geophysics, solar energy, electronic cooling, nuclear 
energy and other industrial system in various sectors. A major 
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limitation against enhancing the heat transfer in such engineering 
systems is the inherently poor thermal conductivity of conventional 
fluids, including oil, water and ethylene glycol mixture. Therefore, for 
more than a century since Maxwell’s theories in 1873, scientists and 
engineers have made great effort to break this fundamental limit by 
dispersing millimeter-or micrometer sized particles in liquids. How- 
ever, the major problem with the use of such large particles is the 
rapid sedimentation of these particles in fluids. 

Maxwell’s concept is old but what is new and innovative in the 
concept of nanofluids is the idea of using nanometer-sized 
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Nomenclature 

Cp specific heat at constant pressure (J kg~! K7!) 
Dr thermal diffusion coefficient (m?/s) 

d diameter (m) 

g gravitational acceleration (m s7?) 

H dimensional cavity height (m) 

k thermal conductivity (W/m K) 

kg Boltzmann constant (J/K) 

p pressure (Pa) 

Pr Prandtl number, Pr = vf /&p 

Pe Peclet number 

Re Reynolds number, Re = PrkeT /37 pF ly 

T dimensional temperature (°C) 

v dimensional nanofluid velocity (m s7!) 


Greek symbols 


u dynamic viscosity (Ns m~?) 


T stress tensor (Pa) 

a thermal diffusivity (m? s7!) 

0) dimensional nanoparticle volumetric fraction 

0) inclination angle 

p density (kg m~?) 

n ratio of the nanolayer thickness to the original 
particle radius 

Subscripts 

f base fluid 

nf nanofluid 

p particle 

k thermal conductivity (W/m K) 


particles to create stable and highly conductive suspensions, 
primarily for suspension stability and for dynamic thermal inter- 
actions. Recognizing an excellent opportunity to apply nanotech- 
nology to thermal engineering, Eastman et al. [1] is the first who 
conceived the novel concept of nanofluid by hypothesizing that it 
is possible to break down these century-old technical barriers by 
exploiting the unique properties of nanoparticles. 

Nanofluids are dilute suspensions of nanoparticles smaller than 
100 nm. The thermal conductivity of the particle materials, metal- 
lic or nonmetallic such as Al,03, CuO, Cu, SiO and TiO, are 
typically order-of-magnitude higher than the conventional fluid. 
Such nanofluids represent a class of its own different from those of 
traditional solid-mixtures, they have superior properties like high 
thermal conductivity [1-5], minimal clogging in flow passage, 
reduced puming power, long term stability and homogeneity. 
Therefore nanofluids have a wide range of potential application 
where improved heat transfer is required. 

The purpose of this paper first is to review the available results of 
the flow and heat transfer presented in literature for natural 
convection of nanofluids, including the experimental and theoretical 
investigations for both single- and two-phase models. The results of 
massive number of different parameters affect the fluid flow and 
heat transfer characteristics are summarized. The second purpose is 
to explore the causes of the enhancement of heat transfer, which 
would be of help in understanding the transport in nanofluids, and 
accurately predict natural convection using nanofluids. 


2. Brief theory 
2.1. Model description 


In the literature, convective heat transfer with nanofluids can 
be modeled using mainly the two-phase or single approach. In the 
two-phase approach, the velocity between the fluid and particles 
might not be zero [6] due to several factors such as gravity, 
friction between the fluid and solid particles, Brownian forces, 
Brownian diffusion, sedimentation and dispersion. In the second 
approach, the nanoparticles can be easily fluidized and therefore, 
one may assume that the motion slip between the phases, if any 
would be considered negligible [7]. The latter approach is simpler 
and more computationally efficient. 


2.2. Conservation equations 


2.2.1. Single phase model 

Although nanofluids are solid-liquid mixtures, the approach 
conventionally used in most studies of natural convection handles 
the nanofluid as a single-phase (homogenous) fluid. In fact, due to 
the extreme size and low concentration of the suspended nanopar- 
ticles, the particles are assumed to move with same velocity as the 
fluid. Also, by considering the local thermal equilibrium, the solid 
particle-liquid mixture may then be approximately considered to 
behave as a conventional single-phase fluid with properties that are 
to be evaluated as functions of those of the constituents. 

The governing equations for a homogenous analysis of natural 
convection are continuity, momentum, and energy equations 
with their density, specific heat, thermal conductivity, and visc- 
osity modified for nanofluid application. The specific governing 
equations for various studied enclosures are not shown here and 
they can be found in different references [18,32,33,55]. 


2.2.2. Two-phase model 

Several authors have tried to establish convective transport 
models for nanofluids. Nanofluid is a two-phase mixture in which 
the solid phase consists of nano-sized particles. In view of the 
nanoscale size of the particles, it may be questionable whether 
the theory of conventional two-phase flow can be applied in 
describing the flow characteristics of nanofluid. On the other hand, 
several factors such as gravity, friction between the fluid and solid 
particles and Brownian forces, the phenomena of Brownian diffusion, 
sedimentation, and dispersion may affect a nanofluid flow. Conse- 
quently, the slip velocity between the fluid and particles cannot be 
neglected for simulating nanofluid flows. Since the two phase 
approach considers the movement between the solid and fluid 
molecule, it may have better prediction in nanofluid study. To fully 
describe and predict the flow and behavior of complex flows, 
different multiphase theories have been proposed and used. The 
large number of published articles concerning multiphase flows 
typically employed the Mixture Theory to predict the behavior of 
nanofluids [8-10]. 

A comprehensive survey of convective transport in nanofluids 
was made by Buongiorno [11], using a model in which Brownian 
motion and thermophoresis are accounted for. Buongiorno 
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developed a two-component four-equation non-homogeneous 
equilibrium model for mass, momentum, and heat transfer in 
nanofluids. The nanofluid is treated as a two-component mixture 
(base fluid+ nanoparticles) with the following assumptions: 


No chemical reactions. 

Negligible external forces. 

Dilute mixture (¢=1). 

Negligible viscous dissipation. 

Negligible radiative heat transfer. 

Nanoparticle and base fluid locally in thermal equilibrium. 


Invoking the above assumptions, the following equations 
represent the mathematical formulation of the non- homogenous 
single phase model for the governing equations as formulated by 
Buongiorno [11]: 


2.2.2.1. Continuity equation. 
V-v=0 (1) 


2.2.2.2. Nanoparticle continuity equation. 
fv-V6=V- (DV +r) (2) 


Here ¢ is nanoparticle volume fraction, Dg is the Brownian 
diffusion coefficient given by the Einstein—-Stokes’s equation: 
kgT 


e= Saud, (3) 


where u is the viscosity of the fluid and dp is the nanoparticle 
diameter and Dr is the thermophoretic diffusion coefficient, 
which is defined as 


H k 
D= (£) (0.26 sexe) 9 (4) 


In Eq. (5), k and kp are the thermal conductivity of the fluid and 
particle materials, respectively. 


2.2.2.3. Momentum equation. 


v-Yv=- | Vp+Y-t+g (5) 
Png 

where 

T=—My(Vu+(Vv)') (6) 


where the superscript ‘t’ indicates the transpose of Vv. 


2.2.2.4, Energy equation. 
PpCp 


v-VT = V(on¢ VT) + 
( nf Pog Cng 


(7) 


(2v4 -YT +D D 


T 


This nanofluid model can be characterized as a ‘two-fluid’ 
(nanoparticles + base fluid), four-equation (2 mass+ 1 momentum + 
1 energy), non-homogeneous (nanoparticle/fluid slip velocity 
allowed) equilibrium (nanoparticle/fluid temperature differences not 
allowed) model. 

Note that the conservation equations are strongly coupled. 
That is, v depends on ¢ via viscosity; ¢ depends on T mostly 
because of thermophoresis; T depends on @ via thermal con- 
ductivity and also via the Brownian and thermophoretic terms in 
the energy equation: ġ and T obviously depends on v because of 
the convection terms in the nanoparticle continuity and energy 
equations, respectively. 


2.3. Physical properties of the nanofluids for single-phase model 


Base nanofluid properties have been published over the past 
few years in literature. However, only recently have some data on 
temperature-dependent properties been provided, even though 
they are only for nanofluid effective thermal conductivity and 
effective absolute viscosity. 


2.3.1. Density 

In the absence of experimental data for nanofluid densities, 
constant-value temperature independent values, based on nano- 
particle volume fraction, are used: 


Par =A—-0) P+ PPp 


2.3.2. Specific heat capacity 

It has been suggested that the effective specific heat can be 
calculated using the following equation as reported in Refs. 
[12-15] as 


(Cp) ng =(1 —$\(Cp )p + P(Cp )p 


Other authors suggest an alternative approach based on heat 
capacity concept [7,16]: 


(PCp)rf = a —)(PCp)p+ P(PCp)p 


These two formulations may of course lead to different results 
for specific heat. Due to the lack of experimental data, both 
formulations are considered equivalent in estimating nanofluid 
specific heat capacity [14]. 


2.3.3. Dynamic viscosity and thermal conductivity 

In this work, different nanofluid models based on a combina- 
tion of the different formulas for the dynamic viscosity and 
thermal conductivity adopted in the studies of natural convection 
are summarized in Tables 1 and 2. 


3. Numerical investigations 
3.1. Single phase model 


Khanafer et al. [17] were the first who studied numerically, the 
buoyancy driven heat transfer enhancement of nanofluid in a 
square cavity. They used Wasp’s model [18] for thermal con- 
ductivity and Brinkman model [19] for viscosity. They illustrated 
that the suspended nanoparticles substantially increase the heat 
transfer rate at any given Grashof number. In addition, they 
indicated that the nanofluid heat transfer rate increases with an 
increase in the nanoparticles volume fraction. They also showed 
that the presence of nanoparticles in the fluid is found to alter the 
structure of the fluid flow. 

Jang and Choi [20] examined the thermal characteristics of free 
convection in rectangular cavity heated from below with nano- 
fluid such as water based nanofluids containing 6 nm copper and 
2 nm diamond, using Jang and Choi’s model [21] for the effective 
thermal conductivity. They concluded that the nanofluids is more 
stable and have larger heat transfer coefficient compared with a 
classical fluid in a rectangular cavity with free convection. 

Kim et al. [22] analyzed the convective instability driven by 
buoyancy and heat transfer characteristics of nanofluids. They 
chose Wang et al. [23] model based on the mean field approach 
for expressing the thermal conductivity enhancement. They 
indicated that as the density and heat capacity of nanoparticles 
increase and the thermal conductivity and the shape factor of 
nanoparticles decrease, the convective motion in a nanofluid sets 
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Table 1 
Various viscosity models for nanofluids. 


Models Dynamic viscosity 
-i Lf 
Brinkman model [19] Lng = ae 
Einstein model [27] Unf = MpA+2.5p) 6<0.05 
Pak and Cho’s Correlation [28] Unf = by +39.1 16+533.9¢7) 


Nguyen et al. model [36] (fitted by Abu-Nada [34]) 


Jang et al. model [41] 


Koo and Kleinstreuer [48] 


Hcuo = —0.6967 + 


0.2597 —30.88 a 


1356.14 
+1.2386+ 7 


15.937 
T 


19652.74 


2 
3, CA Q 
T T +0.01593ġ” +4.38206 T` 147.573 7 


Hano, = €Xp(3.003 —0.04203T—0.5445¢ + 0.0002553T? +0.0524ġ°—1.622¢7') 


Unf = Hp +2.50) 


dp)” 9/3 
1+n(Z) pP (e+1) 


kgT 


=5 x 10bopr |32 -fT, p) 
Ug =5 x 10° Bopp To; Q 


bB=0.0137100ġ) 2829 for ġ<1% 


B=0.0011100ġ) 772? for ġ>1% 


f(T, p) = (—6.04¢ +0.4705)T +(1722.3ġ—134.63) 


1%<h<4 
Maiga model [81] Lng = h(l + 
Brownian model [82] Unf = Hyp 4 
Nguyen model [83] Unf = Mp4 


Gherasim et al. [84] 


% and 300k<T<325k 


+7.3+123¢7) 


+2.5p+6.177) 
+ 0.025 +0.015¢7) 


Unf _ 9. 904¢148¢ 


Uy 


in easily. They also showed that heat transfer coefficient of a 
nanofluid is enhanced by all parameters with respect to the 
volume fraction. 

Jou and Tzeng [24] numerically studied the natural convection 
heat transfer enhancements of nanofluids within a two-dimensional 
enclosure. They analyzed heat transfer performance using Khana- 
fer’s model for various parameters, like volume fraction, Grashof 
number, and aspect ratio of the enclosure. Results showed that 
increasing the buoyancy parameter and volume fraction of nano- 
fluids cause an increase in the average heat transfer coefficient. 

Polidori et al. [25] analytically investigated the natural con- 
vection heat transfer of Newtonian nanofluids over a vertical 
semi-infinite plate. In their study, the dealt with y-Al,03/water 
nanofluids whose Newtonian behavior was experimentally con- 
firmed for particle volume fractions less than 4%. It was shown 
that special care has to be exercised in drawing generalized 
conclusions about the heat transfer enhancement. They found 
that natural convection heat transfer is not solely characterized 
by the nanofluid effective thermal conductivity and that the 
sensitivity to the viscosity model used seems undeniable and 
play a key role in the heat transfer behavior. 

Hwang et al. [26] theorically investigated the thermal char- 
acteristics of natural convection in a rectangular cavity heated 
from below with water-based nanofluids containing alumina 
(Al,03). They used thermal conductivity model of Jang and Choi 
[21] and various models for the effective viscosity. It was shown 
that the experimental results are put between two theoretical 
lines derived from Jang and Choi’s conductivity model used along 
with both Einstein’s model [27] and Pak and Cho’s correlation [28] 
for viscosity. They showed that water-based Al203 nanofluids is 


more stable than base fluid in a rectangular cavity heated from 
below as the volume fraction of nanoparticles increases, the size 
of nanoparticles decreases, or the average temperature of nano- 
fluid increases. They also showed that the ratio of heat transfer 
coefficient of nanofluids to that of base fluid is decreased as the 
size of nanoparticles increases, or the average temperature of 
nanofluids is decreased. In addition, they stated that the effective 
thermal conductivity models used by Kim et al. [22] and Khanafer 
[17] cannot predict the effect of the size of nanoparticles and the 
average temperature of nanofluids on the enhancement of the 
effective thermal conductivity, although the models are able to 
predict the effect of the particle fraction of nanoparticles on the 
enhancement of the effective thermal conductivity of nanofluids. 

Ho et al. [29] aimed to clarify and present a deep explanation 
on the influences of uncertainties due to adopting various 
formulas for the effective thermal conductivity and dynamic 
viscosity. They presented a numerical analysis for natural con- 
vection of water Al2O3 nanofluid in a square vertical enclosure. 
They observed that a significant difference between enhance- 
ments in the dynamic viscosity estimated from the two adopted 
formulas leads to contradictory heat transfer efficacy of the 
nanofluid, so that the heat transfer across the enclosure can be 
found to be enhanced or migtated with respect to the base fluid. 
According to the authors, the contradictory behavior dictated by 
the significant difference between the dynamic viscosity enhance- 
ments calculated, using two formulas might contribute to 
explaining the disparate findings among the numerical predic- 
tions [17] and experimental results [76,77]. 

Oztop and Abu-Nada [30] investigated heat transfer and fluid flow 
due to buoyancy forces in a partially heated enclosure using different 
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Table 2 
Various thermal conductivity models for nanofluids. 


Model 


Wasp model [18] 


Jang and Choi model [21] 


Bruggeman model [23] 


Chon et al. model [35] 


Koo and Kleinstreuer [48] 


Charuyakorn et al. [85] 


Staionary model [86] 


Yu and Choi [87] 


Patel et al. [88] 


Mintsa et al. [89] 


Thermal conductivity 


kaf kp +2kp—2h(kp—kp) 
kp kp+2kf + b(kp—kp) 


3d, 


k 
T = 1) + Bk +18 x 10° TZ kyRej,,, Pr 
f 


nano 


kup — (3@—1)kp/kp +130 —)—1] + VB 


Ag =[(86—-1kp /kp +8 —)—-1) +8kp/kf 


k oreao / dX 93690 fog 07476 
E —1464.7¢° (£) (2) Pr2.9955 Reh 2321 
ky dp kp 


kp + 2k —2p(kp—kp) 


k 
k KEE E 


T 
l E5 x 10 Bboy Co) | FEAT.) 
PPp 


krf kp +2kf—2ġ(kf—kp) 
kp kp+2kf+ġ(kf—kp) 


(1 +boPem) 


ke fi} Moody _ 
ky kj —@)dp 


kup _ Kp + 2kp—2gp(ky—kp)(1 +1? 
kf ky +k + (kp ky) +? 


4 
t 


Kp kpdy 
1 
ky * Rdp—@) 


2kgTdp 
nap updi 


k 
E —1.726+1.0 
f 


nanofluids. The flush mounted heater was located to the left vertical 
wall with a finite length. The temperature of vertical wall was lower 
than that of heater while other walls were insulated. Calculations 
were performed for different Rayleigh numbers, height of heater, 
location of heater, aspect ratio and volume fraction of nanoparticles. 
They found that heat transfer increases with increasing Rayleigh 
number, height of heater and volume fraction of nanoparticles. In 
addition, the heat transfer enhancement was found to be dependent 
on the type of nanofluid and more pronounced at low aspect ratio 
that at high aspect ratio. 

Ögüt [31] studied numerically the natural convection heat 
transfer of water-based nanofluids in an inclined enclosure with a 
constant flux heater. Five types of nanoparticles were taken into 
consideration: Cu, CuO, Al2O0, and TiOz. The results showed that 
the average heat transfer rate increases significantly as Rayleigh 
number and particle volume fraction increase, the effect is more 
pronounced for Ag and Cu nanoparticles. The results also showed 
that the average heat transfer decreases with an increase in the 
length of the heater. As the heater length is increased, the average 
heat transfer rate starts to decrease for a smaller inclination angle. 

Abu-Nada and Oztop [32] analyzed numerically the effect of 
inclination angle on natural convection heat transfer and fluid flow 
in a two dimensional enclosure filled with Cu-nanofluid. It was found 
that the effect of nanoparticle concentration on Nusselt number is 
more pronounced at high volume fraction than at low volume 
fraction. It was also found that the Effects of inclination angle on 
percentage of heat transfer enhancement become insignificant at low 
Rayleigh number but it decreases the enhancement of heat transfer 
with nanofluid. 

Abouali and Falahatpisheh [33] carried out a numerical study 
of internal free convection of Al,O3 water nanofluid in vertical 
annuli. Vertical walls are maintained at constant temperatures 


and horizontal walls are adiabatic. They developed a correlation 
as a function of Nusselt number of base fluid and particle fraction 
which is a monotonically linear decreasing function of particle 
fraction. In addition, they showed that thermal conductivity 
model of Wasp [18] greatly overestimates the Nusselt number 
of experimental data of Putra et al. [76] and the thermal 
conductivity of Jang and Choi [21] agrees well with the experi- 
ment and shows deterioration of natural convection heat transfer. 

Abu-Nada [34] studied the heat transfer enhancement in 
horizontal annuli using variable properties of AlpO3—-water nano- 
fluid. He used Chon et al. [35] expression for conductivity and 
Nguyen et al. [36] experimental data for viscosity. They observed 
that for Ra>10*, the average Nusselt number was reduced by 
increasing the volume fraction of nanoparticles. However, for 
Ra=10°, the average Nusselt number increased by increasing the 
volume fraction of nanoparticles. For Ra=10*, the Nusselt num- 
ber was deteriorated everywhere around the cylinder surface 
especially at high expansion ratio. However, this reduction is only 
restricted to certain regions around the cylinder surface at 
Ra=10?. They observed also that The Nguyen et al. data and 
Brinkman model gives completely different predictions for 
Ra=10* where the difference in prediction of Nusselt number 
reached 30%. However, this difference was less than 10% at 
Ra=10°, 

Aminossadati and Ghasemi [37] presented a study of natural 
convection in a partially heated enclosure from below and filled 
with different types of nanofluids. The effects on the enclosure 
cooling performance of Rayleigh number, solid volume fraction, 
heat source length and location and the type of nanofluid were 
studied. The results indicated that adding nanoparticles into pure 
water improves its cooling performance especially at low Ray- 
leigh numbers. The type of nanoparticles and the length and 


5368 Z. Haddad et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5363-5378 


location of the heat source proved to significantly affect the heat 
source maximum temperature. 

Abu-Nada et al. [38] presented a study of the heat transfer 
enhancement in a differentially heated enclosure using variable 
thermal conductivity and variable viscosity of AlzO3-water and 
CuO-water nanofluids. The results indicated that for an enclosure 
with unity aspect ratio, the average Nusselt number of Al,03- 
water nanofluid at high Rayleigh number was reduced by increas- 
ing the volume fraction of nanoparticles above 5%. However, at 
low Rayleigh number, the average Nusselt number was slightly 
enhanced by increasing the volume fraction of nanoparticles. 
Also, at high Rayleigh numbers, Cuo-water nanofluids manifest 
a continuous decrease in Nusselt number as the volume fraction 
of nanoparticles is increased. However, the Nusselt number was 
not sensitive to the volume fraction at low Rayleigh numbers. It 
was observed that enclosures, having high aspect ratios, experi- 
ence more deterioration in the average Nusselt number when 
compared to enclosures having low aspect ratios. It was also 
observed that the average Nusselt number was more sensitive to 
the viscosity models than to the thermal conductivity models. 

Lin and Violi [39] analyzed the heat transfer and fluid flow of 
natural convection in a cavity filled with Al,03/water nanofluid 
that operates under differentially heated walls. The thermal 
conductivity and dynamic viscosity of the nanofluid are employed 
by Xu’s model [40] and Jang’s et al. model [41], respectively. The 
heat transfer rates were examined for parameter of non-uniform 
nanoparticle size, mean nanoparticle diameter, nanopartticle 
volume fraction, Prandtl! number and Grashof number. They 
illustrated that the heat transfer characteristics of the nanofluid 
can be enhanced as the mean nanoparticle diameter is decreased 
from 250 to 5nm. They attributed these phenomena to the 
dominant effect of the Brownian motion caused by heat convec- 
tion. However, the heat transfer performance of the nanofluid was 
found to be less significant as the dimensionless total thermal 
conductivity of the nanofluid is close to unity due to the increase 
of nanoparticle sizes. This contradictory effect of nanofluids 
was explained to be mainly caused by the effective dynamic 
viscosity. 

Abu-nada and Chamkha [42] focused on the study of the 
natural convection heat transfer characteristics in a differentially 
heated enclosure filled with Cuo-Eg-water nanofluid. Overall, the 
study indicated that the effects of the viscosity models are 
predicted to be more predominant on the behavior of the average 
Nusselt number than the influence of the thermal conductivity 
models. Also, the enclosure aspect ratio is predicted to have 
significant effects on the behavior of the average Nusselt number 
which decreases as the enclosure aspect ratio increases. 

Abu-nada [43] presented a solution of Rayleigh-Bénard convec- 
tion problem of CuO-water nanofluids flowing in a rectangular 
enclosure. The results were presented over a wide range of Rayleigh 
numbers and volume fractions, using variable thermal conductivity 
and variable viscosity models. It has been shown that For Ra > 103, 
the average Nusselt number was reduced by increasing the volume 
fraction of nanoparticles. However, for Ra= 10, the average Nusselt 
number was enhanced by increasing the volume fraction of nano- 
particles. Abu-Nada explained such deterioration in heat transfer 
observed experimentally and showed that the influence of nano- 
particles elucidates two opposing effects on the energy transport: a 
favorable effect driven by the presence of high thermal conductive 
nanoparticles and undesirable effect promoted by high level of 
viscosity experienced at high volume fraction of nanoparticles, 
where at high Rayleigh number the undesirable effect of high 
viscosity becomes dominant, which causes deterioration in heat 
transfer. The results obtained using variable properties of nanofluid 
were compared with those based on constant property simulations 
and it was found that the temperature influence is small compared 


to the influence of high viscosity brought by the presence of high 
concentration of nanoparticles. 

Corcione [44] presented two general correlations for thermal 
conductivity and viscosity of nanofluids using many available 
experimental data. These two correlations are sensitive to types 
of the base fluid and nanoparticle, diameter of nanoparticle and 
the temperature of the nanofluid. Using these new correlations, 
Corcione [45] investigated the heat transfer features of buoyancy 
driven nanofluid inside rectangular enclosures differentially 
heated at the vertical walls. They assumed that nanofluids behave 
more like single-phase fluid rather than a conventional solid- 
liquid mixture, thus, they extended the heat transfer correlations 
available in literature for single phase natural convection in 
enclosure to nanoparticle suspension, by replacing the thermo- 
physical properties appearing in them with the nanofluid effec- 
tive properties calculated at the average temperature. Results 
showed the existence of an optimal particle loading for maximum 
heat transfer. Specifically, for any assigned combination of solid 
and liquid phases, the optimal volume fraction was found to 
increase slightly with decreasing the nanoparticle size, and to 
increase much more remarkably with increasing both the nano- 
fluid average temperature and the slenderness of the enclosure. 

Ghasemi and Aminossadati [46] presented the results of a 
numerical study on the natural convection in a right triangular 
enclosure, with a heat source on its vertical wall and filled with a 
water-Cuo nanofluid. The viscosity and the thermal conductivity 
of the nanofluid were determined based on the models presented 
by Koo and Kleinstreuer [47,48] in which a Brownian motion term 
was added and valid in the range of 1 < ġ < 4%. It was noticed 
that when Brownian motion is neglected, the heat transfer rate 
continuously increases with increasing the solid volume fraction 
at all Rayleigh numbers, whereas, at high Rayleigh numbers, an 
optimum solid volume fraction can be found, which results in the 
maximum heat transfer when Brownian motion is considered. 

Santra et al. [49] simulated the behavior of heat transfer due to 
natural convection in a differentially heated cavity. They calcu- 
lated the shear stresses using Ostwald—de Waele model [50] for 
an incompressible non-Newtonian fluid. It was observed that the 
heat transfer rate decreases with increase in solid volume fraction 
for a particular Rayleigh number. However, it increases with 
Rayleigh number for a particular solid volume fraction. 

Rashmi et al. [51] performed a numerical study of natural 
convection for Al,O03-water nanofluid inside a horizontal cylinder, 
using FLUENT commercial code. Stationary particle model was 
employed for thermal conductivity. The numerical results showed 
decrease in heat transfer with increase in particle volume frac- 
tion. Numerical predictions were also presented to be compared 
with experimental results of Putra et al. [76]. The comparison 
showed a good agreement with the experimental findings. 

Das and Ohal [52] conducted a numerical investigation of 
natural convection inside a square cavity containing Cu—water 
nanofluid with half-active and half-insulated vertical walls. They 
considered five different cases based on the position of the active 
walls. It was found that the nanoparticles when immersed in a 
fluid are capable of increasing the heat transfer capacity of the 
base fluid and this effect is more pronounced for high solid 
volume fraction. 

Ghasemi et al. [53] investigated the effects of magnetic field on 
natural convection in a nanofluid-filled square enclosure for 
different Hartmann number, Rayleigh number and nanoparticle 
volume fraction. They indicated that the heat transfer rate 
increases with an increase of the Rayleigh number but it 
decreases with an increase of the Hartmann number. An increase 
of the solid volume fraction may result in enhancement or 
deterioration of the heat transfer performance depending on the 
value of Hartmann number and Rayleigh numbers. 
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Oztop et al. [54] studied the problem of steady state natural 
convection in an enclosure filled with different types of nano- 
particles (Al,03 and TiOz) numerically. It was observed that the 
addition of nanoparticles into water affects the fluid flow and 
temperature distribution especially for higher Rayleigh numbers. 
In addition, an enhancement in heat transfer rate was registered 
for the whole range of Rayleigh numbers. However, low Rayleigh 
numbers showed more heat transfer increment compared to high 
Rayleigh numbers. It was also noticed that the difference between 
Al,03 and TiOz enhancement becomes more pronounced at high 
volume fractions of nanoparticles. 

Aminossadati and Ghasemi [55] performed a numerical study 
of natural convection in a square cavity filled with a water-CuO 
nanofluid. They considered two pairs of heat source-sink to cover 
the entire length of the bottom wall of the cavity while the other 
walls were thermally insulated. The results showed that regard- 
less of the position of the pairs of source-sink, the heat transfer 
rate increases with an increase of the Rayleigh number and the 
solid volume fraction. 

Jahanshahi et al. [56] examined the influence of uncertainties 
due to adopting various formulas for the effective thermal con- 
ductivity of Silica-water on heat transfer characteristics for 
natural convection in a square enclosure. An experimental setup 
was used to extract the conductivity value of nanofluid. They 
predicted the effective thermal conductivity of the nanofluid using 
the model proposed by Hamilton et al. [57] and experimental 
results. Their findings showed a considerable enhancement in heat 
transfer by increasing the solid volume fraction at any Rayleigh 
number using experimental thermal conductivity, whereas a 
decrease in heat transfer using theoretical thermal conductivity. 

Kefayati et al. [58] reported the numerical results of natural 
convection in enclosures using water/SiOz nanofluid. The main aim 
of their study was to identify the ability of Lattice Boltzmann 
method (LBM) for solving nanofluid and various geometries. The 
simulations were performed for different Rayleigh numbers, solid 
volume fraction and aspect ratio of enclosure. Experimental data 
from Jahanshahi et al. [56] were used to predict thermal conductiv- 
ity of SiO2-water nanofluid. It was concluded that the Nusselt 
number increases with volume fraction for the whole range of 
Rayleigh numbers and aspect ratios. Also the effect of nanoparticles 
on heat transfer augments as the enclosure aspect ratio increases. 

Aminossadati and Ghasemi [59] presented a numerical inves- 
tigation to study the natural convection in an isosceles triangular 
enclosure with a heat source located at its bottom wall and filled 
with an Ethylene Glycol-Cu nanofluid. The results indicated that 
the thermal performance of the enclosure is improved with an 
increase in the Rayleigh number and solid volume fraction. The 
results also indicated that the variation of heat transfer rate with 
respect to the enclosure apex angle, heat source position and 
dimensions is different at low and high Rayleigh numbers. More- 
over, the results obtained from the modified and original Maxwell 
models indicated that the heat transfer rates obtained based on 
the modified Maxwell model are generally higher than those 
obtained based on the original Maxwell model. 

Shahi et al. [60] analyzed three-dimensional flow and heat 
transfer in a single ended tube with non-uniform heat input. They 
assumed that the sealed end of tube to be adiabatic and also the 
tube opening to be subjected to copper-water nanofluid. The results 
indicated that the maximum overall mean Nusselt number is 
obtained at 6=35°, while the maximum output mass flow rate is 
increasing function of the inclination angle. It was also indicated 
that both the Nusselt number and maximum output mass flow rate 
are increasing function of solid concentration, but the presence of 
nanoparticles is more effective at the smaller inclination angles. 

Ghasemi and Aminossadati [61] conducted a numerical study to 
investigate the natural convection cooling of an oscillating heat 


source embedded on the left wall of an enclosure filled with a 
water-based nanofluids (Cu, Al203, TiO2). They examined the effects 
of Rayleigh number, solid volume fraction, heat source position, type 
of nanofluids and oscillation period on the performance of the 
enclosure. The study showed that the addition of nanoparticles into 
the pure water improves its thermal conductivity and enhances the 
heat removal from the heat source. At low Rayleigh number, this 
effect is more noticeable for copper nanoparticles. It was also 
showed that Cu and TiO are the most and the least effective 
nanoparticles, respectively, in the heat removal process. 

Elhajjar et al. [62] used laws of thermodynamics for specific 
heat capacity and thermal expansion coefficient to study numeri- 
cally the Rayleigh-Bénard convection for three types of nano- 
fluids (water-based Al,03, Cu and CuO). It was observed that for a 
fixed value of the fluid Rayleigh number, the nanofluid Rayleigh 
number decreases with the volume fraction of nanoparticles, 
meaning that the nanoparticles delay the onset of convection. 

Sheikhzadeh et al. [63] numerically modeled the buoyancy- 
driven fluid flow and heat transfer in a square cavity with 
partially active side walls filled with Cu-water nanofluid. They 
studied the effects of the position of active portions of the side 
walls, and volume fraction. The result showed that the average 
Nusselt number increases with increasing both the Rayleigh 
number and the volume fraction. Moreover, the maximum aver- 
age Nusselt number for the high and the low Rayleigh numbers 
occur for the bottom-middle and the middle-middle locations of 
the thermally active parts, respectively. 

Saleh et al. [64] studied the heat transfer enhancement 
utilizing water-Cu and water-Al,03 nanofluids in a trapezoidal 
enclosure. Various inclination angles of the sloping wall, solid 
volume fractions, and Grashof numbers were considered and the 
flow and temperature fields as well as the heat transfer was 
analyzed. They found that the structure of the fluid flow within 
the enclosure depends upon, Grashof number, inclination angle of 
sloping wall and nanoparticles concentration and type. In addi- 
tion, the Cu nanoparticles with high volume fraction combines 
with an acute sloping wall are most effective in enhancing 
performance of heat transfer rate. 

Mahmoodi [65] numerically investigated the free convection 
fluid flow and heat transfer of various water base nanofluids. The 
left and right walls of the cavity are maintained at constant 
temperature, while its top and bottom walls are insulated. A thin 
heater is located inside the cavity which its location and length. 
They investigated the effects of the Rayleigh number, position and 
location of the heater, the volume fraction of the nanoparticles, 
and various types of nanofluids on the fluid flow and heat 
transfer. It was observed that for all cases, the rate of heat 
transfer increases with increase in Rayleigh number and the 
volume fraction of the nanoparticles. It was also observed that 
at low Rayleigh numbers the type of nanofluid does not affect the 
heat transfer rate while at high Rayleigh numbers maximum rate 
of heat transfer is obtained by Ag-water nanofluid and its 
minimum is obtained by TiOz—water nanofluid. 

Lai and Yang [66] performed a numerical study to examine the 
natural convective thermal efficacy of Al,03/water nanofluid for 
two circumstances, i.e. at fixed Rayleigh number and at fixed 
temperature difference across the enclosure, using the lattice 
Boltzmann method. Results indicated that the average Nusselt 
number increases with the increase of Rayleigh number and 
nanoparticle volume concentration. In addition, the relative 
Nusselt number increases with an increasing average fluid tem- 
perature and the increase of temperature difference imposed 
between the side walls. The deviations between the predicted 
values of thermal conductivity using various models were found 
not to be remarkable, while significant under-estimations of 
nanofluid viscosity were revealed using different models. 
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3.2. Two phase model 


Two-phase model was proposed to get realistic results. Tzou 
[67,68] analyzed analytically the onset of convection in a horizontal 
layer uniformly heated from below (the Rayleigh-Bénard problem). 
They used a model proposed by Buongiorno [11] for convective 
transport in nanofluid incorporating the effects of Brownian diffusion 
and thermophoresis. In the study of Tzou [67], the case of stress-free 
boundaries was analyzed. In his other work [68], the combination of 
one free and one rigid boundary was investigated, and the case of two 
rigid boundaries was briefly mentioned. On the basis of his calcula- 
tions, Tzou [67,68] concluded that for nanofluids the critical Rayleigh 
number (the value of the Rayleigh number at which instability 
appears) was lower by one or two orders of magnitude than for 
regular fluids. 

The Cheng-Minkowycz problem was solved by Nield and 
Kuznetsov [69] for natural convective boundary-layer flow in a 
porous medium saturated by a nanofluid. They included Brow- 
nian and thermophoresis effects and similarity solution was 
presented. In another study [70], they made an analytical study 
at the onset of convection in a horizontal layer of a porous 
medium saturated by a nanofluid by including Brownian and 
thermophoresis. Again, Nield and Kuznetsov [71] presented a 
work on a linear stability analysis for the onset of natural 
convection in a horizontal nanofluid layer. The similar case is 
studied for natural convective boundary-layer flow of a nanofluid 
past a vertical plate by Kuznetsov and Nield [72]. Khan and Pop 
[73] published a paper on boundary-layer flow of a nanofluid past 
a stretching sheet as a first paper in that field. Their model used 
for the nanofluid incorporates the effects of Brownian motion and 
thermophoresis. They have taken into account the Pr number, 
Lewis number, Brownian motion numbers, N, and thermophor- 
esis number, N;. They indicated that the reduced Nusselt number 
is a decreasing function of each dimensionless number and 
Sherwood number is an increasing function of higher Prandtl 
number and a decreasing function of lower Pr number for each Le, 
N» and N; numbers. Recently, Pakravan and Yaghubi [74] studied 
the thermophoresis, Brownian and Dufour effect on natural 
convective heat transfer of nanofluids simultaneously. They 
showed that the effects of Dufour strongly decreases with 
temperature and Nusselt number of natural convection of nano- 
fluids increases with mean temperature of mixture. 

Recently, Haddad et al. [75] studied numerically the natural 
convection heat transfer and fluid flow of CuO-water nanofluid in 
an open cavity heated from the bottom. A two component non- 
homogenous equilibrioum model was used for the nanofluid, 
incorporating the effects of Brownian motion and thermophoresis. 
Variable thermal conductivity and variable viscosity were taken 
into account. It was showed that higher heat transfer is formed with 
the presence of Brownian and thermophoresis effect. It was also 
found that heat transfer decreases with increasing of nanoparticle 
volume fraction due to increasing of viscosity of the nanofluid. 


4. Experimental investigation 


Experimental works on natural convection in nanofluid enclosure 
are very limited due to difficulties of preparation of nanofluid and 
measurements of the parameters. Putra et al. [76] made an experi- 
mental work on natural convection of AlO; and CuO-water nano- 
fluid inside a horizontal cylinder heated and cooled from the two 
ends respectively. They found that the presence of nanaoparticles 
deteriorates the natural convective heat transfer systematically. They 
observed a systematic degradation of natural convective heat transfer 
with increasing particle concentration characterized by decreasing 
Nusselt number for a given Rayleigh number. The decrease was 


observed to be more severe for aspect ratio of 1.0 than that of 0.5. 
They ascribed the possible reasons of deterioration to the effects of 
particle-fluid slip and sedimentation of nanoparticles. 

Wen and Ding [77] conducted an experimental study of natural 
convective heat transfer of TiO2 (30-40 nm)/water nanofluid in a 
vessel which was composed of two horizontal aluminum disks of 
diameter 240 nm and thickness 10 nm separated by a 10 nm gap. 
The upper surface of the top aluminum disk was exposed to the 
open space and cooled by compressed air. They investigated both 
the transient and steady heat transfer coefficient for different 
concentrations of nanofluids. They found that nanofluids decrease 
the natural convective heat transfer coefficient. Furthermore, such 
deterioration increases with nanoparticle concentrations. They 
attributed the reason of their observation to the particle-particle 
interaction and modifications of the dispersion properties. 

Nnanna [78] presented an experimental system to investigate 
the natural convective heat transfer in a differentially heated 
enclosure, filled with AlpO3—water. It was observed for small volume 
fraction 0.2 < ¢ < 2% the presence of nanoparticles does not impede 
the free convective heat transfer, rather it augments the rate of heat 
transfer. However, for large volume fraction ¢@ > 2%, the convective 
heat transfer coefficient declines due to reduction in the Rayleigh 
number caused by increase in kinematic viscosity. 

Li and Peterson [79] performed an experimental investigation on 
the natural convection heat transfer characteristics of Al203/water 
nanofluids in a vertical cylinder enclosure heated from below. 
Compared with the base fluid (distilled water), the heat transfer 
rate across the enclosure was found increasingly deteriorated with 
the volume fraction of the nanoparticles in the nanofluid. To further 
explain the possible causes, other than the relatively enhanced 
viscosity of the nanofluid, for the deteriorated heat transfer efficacy 
of using the nanofluid, a flow- visualization experiment of sub- 
micron polystyrene-water suspension in rectangular enclosure was 
performed to infer the detrimental influences due to the Brownian 
motions and thermophoresis movements of the particles. 

Ho et al. [80] reported experimental results which illustrated the 
natural convective heat transfer of Al,O3-water nanofluid in a 
vertical square enclosure. The thermophysical properties including 
the particle size, dynamic viscosity, density and thermal conductivity 
were measured as a function of temperature as well as of the 
volumetric particle fraction in the nanofluid. The experimental 
results illustrated systematic heat transfer degradation for the 
nanofluid containing nanoparticles of p > 2% over the entire range 
of the Rayleigh number considered. However, for the nanofluid 
containing much lower particle fraction of 0.1%, a heat transfer 
enhancement of around 18% compared with that of water was found 
to arise in the largest enclosure at sufficiently high Rayleigh number. 

Table 3 shows the summary of published experimental and 
theoretical investigations of the natural convective heat transfer 
performance of various nanofluids. 


5. Recent Studies 


Recently, several studies emerged that analyze natural con- 
vection heat transfer enhancement using nanofluids. Mahmoodi 
[90] proposed a numerical investigation to describe the natural 
convection fluid flow and heat transfer of Cu-water nanofluid 
inside L-shaped cavities. The effects of the Rayleigh number, the 
aspect ratio and the volume fraction were investigated. The 
results indicated that the average Nusselt number, for all ranges 
of cavities aspect ratio, increases with increase in the Rayleigh 
number and the solid volume fraction of the nanofluid. Also, it 
was indicated that the rate of heat transfer was higher for the 
L-shaped cavities with lower aspect ratio. Fattahi [91] investi- 
gated, numerically, natural convection flows in a cavity subject to 
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Table 3 
Summary of experimental and theoretical studies on natural convection of nanofluids. 


Authors Configuration 


Particle Base 
material fluid 


Experimental / 
Numerical 
method 


Parameters 


Observations 


Shahi et al. [15] 


Khanafer et al. 
[17] 


Jang and Choi 
[21] 


L 
Rayleigh- Benard cols Temperature profile 


O@O/ 


PEPE SEED Heating 
Differentially heated cavity 


Kim et al. [22] 


Jou and Tzeng 
[24] 


Polidori et al. 
[25] 


Heated vertical plate 


Hwang et al. [26] 


Ho et al. [29] 


Oztop and Abu- 
nada [30] 


Ogiit [31] 


Abu-nada and 
Oztop [32] 


Cu Water 
Cu Water 
Cu/ Water 
Diamond 

Cu/Ag Water 
Cu Water 


y- AlzO3 Water 


Al203 Water 


Al203 Water 


Al203/ Water 


Cu/TiO2 


Cu/ Water 
Al,03/ 

CuO/Ag/ 

TiO2 


Cu Water 


Finite volume 
method 


Finite-volume 
approach 


Theoretical study 


Theoretical study 


Finite1 difference 
method 


Theoretical study 


Theoretical study 


Finite volume 
method 


Finite volume 
method 


PDQ method 


Finite volume 
method 


10? < qm < 10° 
0<b<5% 


10? <Gr<10° 


0<<25% 
dcy=10 nm 
A=1.0 

300 <T<325K 
dcu=6vnm 
daiamond=2 nm 
p<5% 


0< $< 20% 


10? < Ra < 10° 
0 < p <20% 
dcu=10 nm 
A=0.5, 1.0, 2.0 
0< <4% 

dy- Al203=42 nm 


300 <T <325 K 
10 <d <50 nm 
h <5% 


10? < Ra < 108 
0<ġ<4% 
A=1.0 


10°<Ra<5 x 10° 
0< ġ <20% 
0.5<A<2 


0< ¢ <20% 


104 < Ra < 10° 
0<9<90° 
A=1.0 

0° < p < 120° 
10? < Ra < 107 
0<<10% 
0<@<120° 
A=1 


Both Nusselt number and maximum output 
mass flow rate are increasing function of 
solid concentration. 


The nanofluid heat transfer rate increases 
with an increase in the nanoparticles 
volume fraction. 

Nu = 0.5163(0.4436 + 108?)Gr03123 


The nanofluids are more stable and have 
larger heat transfer coefficient compared 
with a classical fluid. 


Heat transfer coefficient of a nanofluid is 
enhanced by all parameters with respect to 
the volume fraction. 


Increasing the volume fraction of nanofluids 
cause an increase in the average heat 
transfer coefficient. 


The sensitivity to the viscosity model used 
seems undeniable and plays a key role in 
the heat transfer behavior. 

The heat transfer is decreased as the size of 
nanoparticles increases, or the average 
temperature of nanofluids is decreased. 


Significant difference between 
enhancements in the dynamic viscosity 
estimated from the two adopted formulas 
leads to contradictory heat transfer efficacy 
of the nanofluid Nu = C(1 + $)" Ra" 


The heat transfer increases with increasing 
the value of volume fraction of 
nanoparticles. 


The average heat transfer rate increases 
significantly as particle volume fraction 
increases. 


The effect of nanoparticle concentration on 
Nusselt number is more pronounced at high 
volume fraction than at low volume 
fraction. 
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Authors Configuration Particle Base Experimental/ Parameters Observations 
material fluid Numerical 
method 
Abouali and Al203 Water Finite volume 10? < Gr < 10° The thermal conductivity of Jang and Choi 
Falahatpisheh method 0< ¢ġ <6% [22] agrees well with the experiment and 
[33] 1<A<5 shows deterioration of natural convection 
heat transfer. Nuns = Nus(1—cd) 
Abu-nada [34] Al203 Water Finite volume 10? < Ra < 10° For Ra > 104, the average Nusselt number 
method 0< ġ <9% was reduced by increasing the volume 
A=0.2, 0.4, 0.8 fraction of nanoparticles. However, for 
Ra=10°, the average Nusselt number 
increased by increasing the volume fraction 
of nanoparticles. 
Aminossadati Cu/TiO2/ Water Finite volume 10? < Ra < 10° Adding nanoparticles into pure water 
and Ghasemi Al203/Ag approach 0< ġ <20% improves its cooling performance especially 
[37] A=1.0 at low Rayleigh numbers. 
Abu-nada et al. Al20;/ Water Finite volume 10? < Ra < 107 The average Nusselt number was more 
[38] CuO method 0< <9% sensitive to the viscosity models than to the 
1/2<A<2 thermal conductivity models. 
dcuo=29 nm 
daizo3=47 nm 
Lin and violi [39] Al203 Water Weighting 10? <Gr<10° Enhancement and mitigation heat transfer 
function method 5 <dapo3 < 250 nm can be attributed to the dominant effect of 
0<f<5% the Brownian motion caused by heat 
A=1.0 convection. 
Abu-nada and CuO Water/ Finite volume 10? < Ra < 10” The effects of the viscosity models are more 
Chamkha [42] LLLLLLbL Legs EG method 0< ¢ġ <6% predominant on the behavior of the average 
1/2 <A <2.0 Nusselt number than the influence of the 
thermal conductivity models. 
Abu-nada [43] e e CuO Water Finite volume 10? < Ra < 10° For Ra > 103, the average Nusselt number 
method 0<p<9% was reduced by increasing the volume 
ia Š m A=2.0 fraction of nanoparticles. However, for 
ni G por Ra=10?, the average Nusselt number was 
ty ' enhanced by increasing the volume fraction 
[a asic of nanoparticles 
hs 
Corcione [45] Differentially heated enclosure Cu/TiO2/ Water/ Theoretical study 249<T<324K The optimal volume fraction is found to 
Al203 EG 25<d<150nm increase slightly with decreasing the 
1<A<40 nanoparticle size, and to increase much 
p<6% more remarkably with increasing both the 
nanofluid average temperature and the 
slenderness of the enclosure. 
Ghasemi and CuO Water Finite volume 10? <Ra<10° When Brownian motion is neglected, the 


Aminossadati 
[46] 


method 


1<p<4% 
1/2<A<2.0 


heat transfer rate continuously increases 
with increasing the solid volume fraction at 
all Rayleigh numbers. 
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Table 3 (continued ) 


Authors Configuration Particle Base Experimental / Parameters Observations 
material fluid Numerical 
method 
Santra et al. [49] k Cu Water Finite volume 10f<Ra<10 The heat transfer rate decreases with 
method 0<ġ<5% increase in solid volume fraction for a 
dcu=100 nm particular Rayleigh number. However, it 
A=1.0 increases with Rayleigh number for a 
i | 9 particular solid volume fraction. 
Ty Te 
1 
x 
Rashmi et al. [51] y Al203 Water Fluent software œ =1%,4% The rate of heat transfer decreases with 
A=1.0 increase in particle volume fraction. 


Das and Ohal ER E erent Cu Water Finite volume 10f<Gr< 10” The rate of heat transfer increases with an 
[52] method 0< ġ <20% increase in the nanoparticles volume 
| A=1.0 fraction. 
Te T 
? ASSAF 
adiabatic 
Ghasemi et al. > hale Al203 Water Finite volume 10° <Ra <10” The heat transfer rate increases with an 
[53] method 10? < Ha < 107 increase of the Rayleigh number but it 
ERTA 0< <6% decreases with an increase of the Hartmann 
Nanofluid A=1.0 number. 
Oztop et al. [54] WA Cu/TiO2/ Water Finite volume 10° < Ra < 10° An enhancement in heat transfer rate was 
Al203 method 0< <10% registered for the whole range of Rayleigh 
Ao 0<@<90° numbers. 
j Q A=1.0 
y 
N 
$s 
E30- 
33 
he 
Aminossadati CuO Water Finite volume 10° < Gr < 10° The heat transfer rate increases with an 
and Ghasemi case (1) method 0<p<4% increase of the Rayleigh number and the 
[55] A=1.0 solid volume fraction. 
Water-CuO 
Jahanshahi et al. SiO2 Water Finite volume 107 < Ra < 107 The mean Nusselt number increases with 
[56] method/ 0< <4% volume fraction for the whole range of 
Experimental A=1.0 Rayleigh numbers. Although by using the 
dsioz=12 nm theoretical formulations for conductivity no 


enhancement was observed. 


SiOz Water Lattice 10? <Ra<10° The Nusselt number increases with volume 
Boltzmann 0< p <4% fraction for the whole range of Rayleigh 
method 1/2 <A < 2.0 numbers and aspect ratios. 


Kefayati et al. 
[58] 
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Authors Configuration Particle Base Experimental/ Parameters Observations 
material fluid Numerical 
method 
Aminossadati Cu EG Finite volume 10? <Ra<10° The thermal performance of the enclosure is 
and Ghasemi method 0<p<5% improved with an increase in the Rayleigh 
[59] number and solid volume fraction. 
Ghasemi and Cu/TiO2/ Water Finite volume 10? < Ra < 10° The addition of nanoparticles into the pure 
Aminossadati 2 Al203 method 0< ġ <20% water improves its thermal conductivity 
[61] H A=1.0 and enhances the heat removal from the 
z heat source. 
L 
Elhajjar et al. [62] Rayleigh-Bénard configuration | Cu/TiO2/ Water 0<Ra<10* The nanofluid Rayleigh number decreases 
Al203 0< <8% with the volume fraction of nanoparticles. 
Sheikhzadeh Cu Water Finite volume 10? < Ra < 10° The average Nusselt number increases with 
et al. [63] method increasing both the Rayleigh number and 
the volume fraction. 
0<¢ <15% Nuavg = (0.4778? +0.3554)Ra?25 
Saleh et al. [64] Cu/Al,03 Water Finite difference 10° < Gr < 10° The Cu nanoparticles with high volume 
method fraction combines with an acute sloping 
wall are most effective in enhancing 
performance of heat transfer rate. e 
0< <26% k So 
0° <p x26" Nutayg = 0.231Gr/032/ (e) 
Mahmoodi [65] Cu/TiO2/ Water Finite volume 10? < Ra < 10° Acute sloping wall and Cu nanoparticles 
Ag method 0< ġ <9% with high concentration are effective to 
A=1.0 enhance the rate of heat transfer. 

Lai and Yang [66] Al203 Water Lattice 10? < Ra < 10° The average Nusselt number increases with 
Boltzmann 0< p <4% the increase of Rayleigh number Ra and 
method 0.006 < AT < 10.5k nanoparticle volume concentration. 

A=1.0 
Z. Haddad et al. Calaca 7 CuO Water Finite volume 10? < Ra < 10° Higher heat transfer is formed with the 
[75] method 0< p <9% presence of Brownian and thermophoresis 
w i n A=2.0 effect. 
Perea) 8 Pert 
poo 
| 1 Hot Pise (Ta) 
Putra et al. [76] CuO/ Water Experimental 10’ <Ra<10° A systematic degradation of natural 
Al,03 study h = 1%,4% convective heat transfer with increasing 
daiz03= 131.2 nm particle concentration for a given Rayleigh 
dcuo=87.3 nm number. 
A=0.5, 1.0 


1. Cylindrical block 
2, End cover as heating surface, 


6 The piston shaft 
7. Cooling water inlet and outlet 
3. End cover as cooling surface & Narrow tube 

Thermec 


4 


Cap, 9. 
$. Resistance heating elements 
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Authors Configuration Particle Base Experimental / Parameters Observations 
material fluid Numerical 
method 
Wen and Ding TiOz Water Experimental dsio2 = 30-40 nm Nanofluids decrease the natural convective 
[77] study A=0.5, 1.0 heat transfer coefficient, such deterioration 
0O<ġp<1% increases with nanoparticle concentrations. 
10f<Ra< 10° 
Experimental setup: (A) filling vessel, (B) discharging vessel, 
(C) top aluminium plate, (D) bottom aluminium plate, (E) 
insulation material, (F) insulation cover, (G) level indicator, 
(H) heating clement, (I) heat flux sensor, (J) DC power 
supply, (K) data acquisition system (DAQ), (L) computer. 
Nnanna [78] Al203 Water Experimental 1x10’ <Ra<3x 10’ Nusselt increased with small volume 
study fraction,0.2 < @ < 2%, and no significant 
change in Nusselt number is observed in the 
range 2% <p < 7.9%. 
0< ġ <8% Nu = e-Rarpee™) 
daiz03=27 nm 
A=1.0 
Experimental setup. 1) insulation, 2) resistive electrical heater, 
3) heated wall, 4) nanofluid medium, 5) thermocouple wire, 
6) the black dots indicate the thermocouple locations, 7) 
H3S mm eS m 
Li and Peterson Al203 Water Experimental 0.5% < h < 6% The heat transfer rate was found 
[79] study dapo3=47 nm increasingly deteriorated with the volume 
> 8 x 10° < Ra < 2.8 x 104 fraction of the nanoparticles. 
— 
Experimental test facility: (Left) 1: horizontal plane table; 2: 
base plate; 3: vertical support; 4: load screw; $: heat 
dissipater; 6: Thermal insulation wrapper; 7: upper copper 
bar; 8: sample charge tubes; 9: o-ring; 10: lower copper bar; 
11: heat sink; 12: coolant inlet and outlet; (Right) the 
thermocouple distribution on both copper 

Ho et al. [80] Al203 Water Experimental 6.21 x 10° < Ra < 2.56 x 108 Systematic heat transfer degradation for the 

study 0.1% < h < 4% nanofluid containing nanoparticles of 
o = 2% over the entire range of Rayleigh 
number. 

Mahmoodi [90] Cu Water Finite volume 10? < Ra < 10° The average Nusselt number increases with 
method Simpler 0< ¢ġ <10% increase in the solid volume fraction of the 
algorithm A=0.2,0.4,0.6 nanofluid. 

-yy 

Fattahi [91] Differentially heated cavity Cu/Al203 Water Lattice 10? < Ra < 10° The average Nusselt number increases for 
Boltzmann 0<ġ<5% both nanofluids. Further, the effects of solid 
method volume fraction for Cu are stronger than 

Al203. 
Kefayati et al. adiabatic w, Cu Water Lattice 10f<Ra < 10° The average Nusselt number increases with 
[92] ©. > Boltzmann 0< <5% augmentation of Rayleigh number and the 
method 1/2 <A < 2.0 volume fraction of nanoparticles for whole 
3a ranges of aspect ratios. 
Tu $ Te 
aA i 
t adiabatic N) 
Abouali and Six different enclosure shapes CuO/ Water Fluent software - The numerical simulations for natural 
Ahmadi [93] Al203 convection of nanofluids could be obtained 
theoretically in a much simpler way. 

Guo et al. [94] Differentially heated cavity Al203 Water Lattice 10° < Ra < 10° The critical Rayleigh number is lower for 
Boltzmann 0<¢ġ<15% nanofluid than that for base fluid. 
method A=1.0 Nu =0.05Ra?45 


5375 


5376 


Table 3 (continued ) 


Z. Haddad et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5363-5378 


Parameters 


Observations 


Authors Configuration Particle Base Experimental / 
material fluid Numerical 
method 
Mahmoodi and Cu Water Finite volume 
Hashemi [95] method Simpler 
algorithm 
Soleimani et al. Cu Water Finite Element 


Method 


10? <Ra < 10° 
0<p<10% 
A=0.2,0.4,0.6,0.8 


10? <Ra<10° 


The mean Nusselt number increased with 
increase in volume fraction regardless 
aspect ratio of the enclosure. 


The effect of nanoparticles is more 


[96] 


0< <6% pronounced at low Rayleigh number than at 
45°<m<180° high Rayleigh number. 
A=1.0 


different side wall temperatures. The fluid in the cavity was a 
water-based nanofluid containing Al2O3 and Cu nanoparticles. 
The effective thermal conductivity and viscosity of nanofluid were 
calculated using Chon and Brinkman models. The results indi- 
cated that by increasing solid volume fraction, the average 
Nusselt number increased for both nanofluids. Kefayati et al. 
[92] numerically studied natural convection of Cu-water nano- 
fluid in an open enclosure using Lattice Boltzmann method (LBM). 
The results indicated that the average Nusselt number increased 
with augmentation of the Rayleigh number and the volume 
fraction of nanoparticles for the whole range of aspect ratios. 
However, the average Nusselt number decreased as the aspect 
ratio increased. Abouali and Ahmadi [93] analyzed natural con- 
vection of nanofluids in horizontal and titled square, horizontal 
and vertical annuli, a triangular enclosure and the Rayleigh- 
Bénard convection configuration. The study was performed for 
Al203 and CuO nanofluids. The results were compared with the 
prediction of a theoretical approach which uses the Nusselt 
correlations available in literature for pure fluids. Guo et al. [94] 
conducted a numerical study for natural convection heat transfer 
in a square cavity using the LBM. They derived a correlation for 
the mean Nusselt number. Mahmoodi and Hashemi [95] pre- 
sented a numerical study of natural convection fluid flow and 
heat transfer inside C-shaped enclosures filled with Cu-water 
nanofluid using finite volume method. It was found that the mean 
Nusselt number increases with increase in the Rayleigh number 
and volume fraction of Cu nanoparticles regardless of aspect ratio 
of the enclosure. Also, the effect of nanoparticles on enhancement 
of heat transfer at low Rayleigh numbers is more significant than 
that at high Rayleigh numbers. Moreover, it was shown that the 
rate of heat transfer increases as the C-shaped enclosure becomes 
narrower. Soleimani et al. [96] proposed the control volume based 
finite element method (CVFEM) to investigate natural convection 
heat transfer in a semi-annulus enclosure filled with Cu-water 
nanofluid. The inner and outer semi circular walls are maintained 
at constant temperatures while the two other walls are thermally 
insulated. The numerical investigation is carried out for different 
governing parameters: namely the Rayleigh number, nanoparticle 
volume fraction and the angle of turn for the enclosure. The 
effective thermal conductivity and viscosity of nanofluid were 
calculated using the Maxwell-Garnetts (MG) and Brinkman mod- 
els, respectively. It was noticed that the streamlines and iso- 
therms strongly depend on the volume fraction of nanoparticle. 


6. Conclusion 


This paper presented a comprehensive review for the recent 
work published on heat transfer enhancement in natural convec- 
tion using nanofluids. The paper presented experimental, numer- 
ical and analytical published work in literature. The numerical 
work included single-phase and two-phase models. Various 
thermal conductivity and viscosity models employed by various 
researchers were presented. 

The literature survey shows that natural convection heat 
transfer enhancement using nanofluid is still controversial and 
there is an on-going debate on the role of nanoparticles on the 
heat transfer enhancement since there is a lack of extensive 
studies on this topic. Clearly, most of numerical results showed 
that nanofluids significantly improve the heat transfer capability 
of conventional heat transfer fluid. Whereas, experimental results 
showed that presence of nanoparticles deteriorates heat transfer 
systematically. It was observed that the experimental results 
were mostly obtained using one type of nanofluid (Al203-water). 
Thus, benchmark experiments are highly desired to gage the 
validity of the numerical results. Concerning the numerical 
studies which try to explain the observed anomalous enhance- 
ment in heat transfer, it will be necessary to consider not only one 
single phase model but two-phase approach, which seems a 
better model to describe the nanofluid flow since slip velocity 
between the particle and base fluid plays important role on the 
heat transfer performance of nanofluids. 
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